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Abstract 

The finite conductivity corrections to the Casimir force in two configurations are 
calculated in the third and fourth orders in relative penetration depth of electromagnetic 
zero oscillations into the metal. The obtained analytical perturbation results are 
compared with recent computations. Applications to the modern experiments are 
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1. Introduction 

In 1998 just 50 years have passed after the publication of the famous 
paper by Casimir [1] (see the Proceedings [2] especially devoted to this event). 
Accidentally the experimental interest in the Casimir effect was rekindled at 
the same time. In Ref.[3] the torsion pendulum was used to measure the 
Casimir force between Cu plus Au coated quartz plate, and a spherical lens 
in a distance range from 0.6 to 6 fim. The accuracy of order 5% was 
claimed in [3] for the agreement of the measurement with theory. In Refs.[4- 
6] the Casimir force between Al plus Au/Pd coated disk and a sphere was 
measured for surface separations between 0.1 to 0.9 using the Atomic 
Force Microscope. The deviation between theory and experiment was shown 
to be of around 1% at the smallest surface separation [4,6,7]. The obtained 
experimental results and the extent of their agreement with theory were used 
to establish stronger constraints for the parameters of hypothetical long- 
range interactions predicted by the unified gauge theories, supersymmetry 
and supergravity [8-11]. 

To be confident that data fit theory at a level of about several percent, the 
different corrections to the ideal expression for the Casimir force should be 
taken into account. The main contribution is given by the corrections due to 
finite conductivity of the boundary metal, its roughness and due to non-zero 
temperature (see [12] for review). Experimental data of [3] do not support 
the presence of any of these corrections although they should contribute at a 
level of 5%. By contrast, in [4,6] the surface roughness and finite conductivity 
corrections are of great concern (the temperature corrections are negligible 
in the measurement range of [4,6]). 

The subject of the present paper is the calculation of higher order finite 
conductivity corrections to the Casimir force in relative penetration depth of 
electromagnetic zero oscillations into the metal. We consider configurations 
of two plane parallel plates and a sphere above a plate. The first order 
finite conductivity correction was found in [13] for configuration of two plane 
parallel plates with an error in numerical coefficient corrected in [14]. Later 
the correct result was reobtained in [15]. Second order correction was firstly 
found in [16] (see also [12]). It was modified for the configuration of a sphere 
above a disk in [17] by the use of Proximity Force Theorem (PFT) [18]. 
The results of [16,17] for the Casimir force up to the second power in relative 
penetration depth are in common use when discussing the recent experiments 
(see, e.g., [4,6,7,19-21]). In [7,10] the third and the fourth order corrections 
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were obtained approximately from the interpolation formula. They allowed 
to achieve the excellent agreement between theory and experiment. 

In [19] numerical calculation of the Casimir force with account of finite 
conductivity has been attempted based on the tabulated data for the comp- 
lex dielectric permittivity as a function of frequency. The same computation 
was repeated in [21] with the diverged results. The reason of these differences 
was interpreted in [21] as the invalid manipulation of optical data in [19]. Our 
analytical calculation of higher order conductivity corrections agrees with the 
results of [21] in the application range of perturbation approach. As shown 
below the perturbation results obtained in the context of plasma model are 
valid with rather high accuracy when the distance between the test bodies 
is larger than the plasma wavelength (not much larger as advocated in [20]). 
This gives the possibility, in some instance, to use the plasma model for the 
distances of order or even less than the characteristic absorption wavelength 
of test body material. 

The paper is organized as follows. In Sec. 2 the general finite results for 
the Casimir energy density and force are briefly presented. Sec. 3 contains 
derivation of the third and fourth order conductivity corrections. In Sec. 4 
the obtained perturbation results are compared with numerical calculations. 
Sec. 5 contains conclusions and discussion. 



2. Casimir energy density and force between realistic 
materials 

Let us consider two semi-infinite solids with dielectric permittivity e(u) 
separated by a plane-parallel gap of width a. The surfaces of the bodies 
are planes z = 0, a. The Casimir energy density and force acting between 
these bodies can be found most simply following [22,23] (see also [24] for the 
multilayered walls). Some additional clarification is given here in the case in 
which the finite energy density rather than force is the subject of interest. 

The mode frequencies of electromagnetic field are found from Maxwell 
equations supplemented by the standard boundary conditions at z = 0, a. 
Two types of such frequencies u£ ,2) (surface modes [25]) corresponding to two 
polarizations of the electric field are the solutions of the equations 



A«(w£ } r ) = e 2 K^e- K '- a \{K e + eK) 2 e Ka - {K e - eK) 2 e 



Kt 



= 0, 



ti 2) {^\) = K 2 e e- K * a \{K e + Kf e Ka - (K e - K) 2 e~ Ka ] =0. (1) 
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Here the following notations are introduced 



K 2 = k 2 - K 2 = k 2 - e^, 



u 



(2) 



and k is the two-dimensional propagation vector in the x?/-plane. 

The infinite zero-point energy of the electromagnetic field between the 
plates is given by [23,24] 



fc,n 



,( 2 ) 



(3) 



where the sum in continuous index k is actually an integral. Introducing the 
length L for the x, y sides of the plates we obtain the vacuum energy density 



(2) 
n I i 



(4) 



which is also infinite. 

Summation in (g) over the solutions of ([!]) can be performed with the help 
of the argument principle 



2X 



(1,2) 



1 



2m 



J cudlnA^ 2 \cj) + JtudlnA^ 



(5) 



where contour C + is a right semicircle of infinite radius in complex cj-plane 
with a center at the origin (note that the functions A^' 2 ^(o;) have no poles). 
The integral over C + can be simply calculated when it is considered that 



lim e(u) = 1, 



lim — ^ = 



(6) 



along the arbitrary radial direction in complex plane. The result (infinite) 
does not depend on a. It is given by 



JtudlnA^ 2) {u) = 3 J dcu. 



(7) 



c+ 



Introducing a new variable £ = —icu in (§), (|7j) and performing a partial 
integration one obtains 



1 00 ^ 

n -oo Ci 



(8) 
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Now let us turn to the removing of divergencies (this important point was 
not discussed in [22-24]). It is apparent that for the infinitely remote plates 
the regularized physical vacuum energy density should vanish [12]. In the 
limit a — > oo we have from (181) 



\n 'Joe 27r_4 27r d 

where by the use of ([I]), it follows 

= e 2 K 2 £ e- K ^ a (K £ + eK) 2 e Ka , = K 2 e eT K ^ (K £ + K) 2 e Ka . 

For a regularized quantity the result is 



E 



(1.2) A _ 



reg 



-oo 



A (1 ' 2) (^) 



(9) 



(10) 



(11) 



where we are guided by the argument of infinitely remote plates. 

Substituting the regularized quantities fiTT]) into (g) instead of (§) we 
obtain the final expression for the Casimir energy density between plates 



oo oo 

£reg(a) = — 2 / kdk J d£ [in A' 1 ' (if) + In A< 2 >('<f )] 



(12) 



where 



A(i) = 



AW 

oo 



= 1 



oo 



(Ks 



(K e + Kf 



K)\-2Ka 

' 1 



(13) 

and also use was made of the fact that A^ 1 ' 2 ) are even functions of £. 

Notice that in [23] no finite expression for the energy density was obtained. 
In [24] the omission of infinities was performed implicitly without a physical 
justification. To illustrate this, in [24] instead of Eqs.fll]) the result of their 
division by the terms containing exp(i^a) was used in spite of the fact that 
on C + such operation is the division by infinity. Fortunately, this operation 
did not influence the final result for the energy density obtained in [24] which 
is perfectly correct. 

One can obtain the Casimir force between plates from ( |T2|) 

Or> f \ j- OO OO ( V / -rr T^\0 1 — 1 

_ dS reg {a) 





kdk Jd^K 



Ks + sKf e2Ka 



+ 



{Ke±Kf 2Ka 

(K £ - Kf 



i 



1 



1 



(14) 
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which is exactly equivalent to Lifshitz result [23,26]. Using PFT [18] we 
obtain from (\T2\ ) the Casimir force acting between a plane plate and a spherical 
lens or a sphere which is given by 



oo oo 



Fi(a) = 2irRS reg (a) = ^JkdkJd£ [in A (1) (^) + In A (2) (^)] • (15) 

o o 



Both Eqs.(P^) and (|T5|) are used below to calculate higher order conductivity 
corrections to the Casimir force between realistic metals in two configurations 
under consideration. 



3. Higher order conductivity corrections 

It is common knowledge that the dominant contribution to the Casimir 
force comes from frequencies £ ~ c/a. We consider the micrometre domain 
with a from a few tenths of a micrometre to around a hundred micrometers. 
Here the dominant frequencies are of visible light and infrared optics. In 
this domain, the plasma model works well and the dielectric permittivity of 
a metal can be presented as 



e(uj, 



1 



UJ 



eft) 



l + > 



(16) 



where the plasma frequency cu p is different for different metals. 

The case of plane parallel plates will be our initial concern. Introducing 
new variables p and x according to 



k 2 = w 



1 



Cr 



cx 



2pa 



(17) 



we transform Eq.([T3j) into the form 




(s + pef 
-pe)' 



-i 



+ 



s + p)' 
(s-p)' 



where 



(18) 

s = ye - 1 +p 2 - (19) 

Let us expand the expression under the integral with respect to p in powers 
of a small parameter 

£ c x 5q x 
p a 2p' 



a — 



2u p a 



(20) 
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where 5q = X p / (2tt) is the effective penetration depth of the electromagnetic 
oscillations into the metal. Note that in terms of this parameter e{uS) = 
l + (l/a 2 ). 

After the straightforward calculations one obtains 

■l 



v s +pey 
[s — pe)' 



1 



2A 



1 



1 



1 



4A 
P 

1 2 , n J 



8A 



a + — (2A 



lW 



P 



+ —r(-6 + 32A - 32A Z + 2p z - p 4 W 
p 6 

8A,„ , jW „ _ , _ „ 2 n J , ,.4\ 4 



:2i: 



-(2,4 - 1)(2 - 16A + 16,4' - 2/ + j9 4 )a 4 + O^) 



where A — e x / (e x — 1). 

In perfect analogy, the other contribution from (|TB1 ) is 

-l 



[s + pY 
[s-p)'- 



1 



1 



1 



2_2 



1 - 4Apa + 8,4(2,4 - l)p l a 



+ 2A(-h + 32,4 - 32,4 2 )p 3 a 3 

+8,4(1 + 18,4 - 48,4 2 + 32^ 3 y a 4 + 0(a 5 ; 



(22) 



(note that this expression actually does not depend on p due to (|20|)). 

After substitution of (ggj ) into (|T5|) all integrals with respest to p have 
the form Jq 00 dpp~ k with k >2 and are calculated immediately. The integrals 
with respect to x have the form 



00 



(e 1 - l) m+1 



(23) 



and can be easily calculated with the help of [27]. Substituting their values 
into ( |I5|) we obtain after some transformations the Casimir force between 
metallic plates with finite conductivity corrections up to the fourth power in 
relative penetration depth 



F v {d) = Ff\a) 



1 



640 



3 a 
2800 



9 



1 



a 

1637T 



7 

2\ x4" 
°0 



1 



7T 



210/ a 3 



7350 



a 4 



(24) 



where Fj )(a) = -(A) /(240a 4 ). 

As was mentioned in the Introduction, the first order correction in ([24]) 
was obtained in [13-15]. The second order correction was obtained in [16] 



8 



(see also [12]). The third and fourth order corrections which are obtained 
here are important for the recent Casimir force measurements (see Sec. 4). 

Now let us turn to the configuration of a lens or a sphere above a plate. 
Introducing the new variable ( [T7|) into (|Tol) we get the Casimir force 



F,(a) 



TicR 7 2 i °?dp 

x ax / — r < In 



167ra 3 



p* 



pe)' 



[s + pe)' 



+ In 



1 



[s + p) : 



(25) 

Bearing in mind the further expansions it is convenient to perform in (|Z5p 
integration by parts with respect to x. The result is 



FAa) = 



hcR 
48ira 3 



OO OO 7 

1 r 



[s + pe) 2 e x — (s — pe) 2 



(s + p) 2 e x — (s — p) 2 



(26) 



The expansion of the first term under the integral in powers of the parameter 
a introduced in ( |2D| ) is 



(s + pe) 2 e x — (s — pe) 2 



1 f 4 

1 H (1 - Ax)a 



8-4 , „ nA s2 2 

+ -5-(-2 - x + 2Ax)oT + -r- 
p z x p 6 x 



1 L px 

2 - 6p 2 + 3/ 



+Ac(-6 + 32,4 - 32,4 2 + 2p 2 - p A ) + 16,4(2,4 - 1)] a 3 
8A 

+ -- \-8 + 32,4 - 32,4 2 + 8/ - 4p 4 

+x(2,4 - 1)(2 - 16,4 + 16,4 2 - 2/ + j9 4 )l a 4 + 0(a 



(27) 



In the same way for the second term under the integral of flZEl ) one obtains 



(s - p) 2 - (s + p) 



2 d [s-pf 
dx (s+p) 2 



1 



1 



[s + p) 2 e x — (s — p) 2 e a 
+ —(-2 - x + 2Ax)p 2 o! 2 + - 



l + -(l-AcW* (28) 
- 1 - 16,4 + 32,4 2 - 5Ax 



x 



x 



+Z2A l x - ?>2A 6 x) p 3 a 3 + — ( -4 + 32,4 - 32,4 



8,4 



x 



x 



+18Ax - 48A 2 x + 32A 3 x) p A a A + 0(a 5 ) 
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Substituting (|27| ), ( p8|) into ( |25Q we firstly calculate integrals with respect 
to p. All integrals with respect to x are of the form (^3|). Calculating them 
we come to the following result after long but straightforward calculations 



Fi(a) = F/ 0) 



a) 



1 



4 ^o + M 
a 5 a 2 

400 




7T 



210 



a 3 



(29) 



where (a) = — (7r 3 foci?)/(360a 3 ). Note that the first order correction from 
(|29|) was firstly published in [3] and the second order one in [17]. 

Although the results ([M]) and ( |2D1 ) for two configurations were obtained 
independently they can be tied by the use of PFT. By way of example, the 
energy density associated with the fourth order contribution in ([M]) is 



4 4) (a) = / F^\a)da 



57T 2 flC 



1 



163tt 2 \ S% 



a 



T 



(30) 



27 \ 7350 

Then the fourth order contribution to the force between a plate and a lens 



F l [4 \a) = 27rRE^{a) 



27a 3 



1 



163tt 2 \ 8$ 
7350 j a 1 



(31) 



agrees with 
way. 



The other coefficients of ( |2DQ can be verified in the same 



4. Comparision with numerical calculations 

In this section we consider the application range of the expressions 
and (|29| ) for the Casimir force which take into account higher order conductivity 
corrections. It is apparent that the greater the distance a between the test 
bodies, the more exact are the perturbation formulas obtained up to the 
fourth power in small parameter 5o/a. Let us compare the correction to 
the Casimir force between two plane parallel plates given by Eq. (^§J) with 
the numerical results. These results were obtained in |JTj for three metals 



(Au, Cu and Al) by the numerical integration of the formulas which are 
equivalent to (|T%| ). In doing so the tabulated data [28] for the complex 
dielectric permittivity was used. The quantity e(i£) was obtained through 
the imaginary part of dielectric permittivity by the use of dispersion relation 

m]. 
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In Fig. la, the solid line represents computational results of Ref. Q2Tj] for 
Fp/F^ in case of Al depending on distance between the plates a. The short- 
dashed line is obtained from Eq. with the value of plasma wavelength 
X Al = 98 nm; the long-dashed line takes account the terms of ( |Z31 ) up to 
the second power only. It is seen that (121 ) is in excellent agreement with 
computational results of [0T| for all a > X Al . For example, for a = 0.1 //m, 
0.5 /im and 3 /zm it follows from (|Zj) that F p /F^ = 0.56, 0.85 and 0.97 which 
can be compared with computations of [ETf: 0.55, 0.85, and 0.96, respectively. 

In Fig. lb, the analogical results for Cu and Au are shown. The dashed 
lines were obtained with X^ Au = 132 nm. For the typical distances indicated 
above it follows from that F p /F^ = 0.60, 0.81 and 0.96 which can be 



compared with the values: 0.48, 0.81 and 0.96 [jZT]l . The difference in the first 



values is due to Xp U,Au > 100 nm, i.e. (^§J) is not applicable for a = 100 nm 
in case of Cu and Au. For a > Xp U,Au the results agree perfectly well. Note 
that the values of plasma wavelength X p = c^Ttm/ (eyN), where m is the 
effective mass of conduction electrons, N is their density are known not very 
precisely. For Al, usually Xp = 100 nm is used |28j| . For Au and Cu the 



value X^ Au = 136 nm was estimated recently [|2Ifl . We used a bit different 



values which provide the smallest rms deviation between the computational 
results and the ones obtained from (|2j) (in [|29j X^ u = 132 nm). The values 
of Fp/F^ at typical distances do not depend on the change of X p for 2-3 
percent which is the uncertainty of the current information regarding X p . 

Now let us turn to the Casimir force between a plate and a lens. The 
numerical results were obtained in [|21j by the integration of equation equivalent 
to ([Z5]). In Fig. 2a, the results for Al bodies are shown, and in Fig. 2b - 
for Cu or Au ones. Solid lines represent computations of [|2T|, short- and 
long-dashed ones are obtained from Eq. ( |2T?| ) used in full or up to the second 
power terms. In both figures the fourth-order perturbation results are in 
excellent agreement with computations for all a > X p . At the distances 
a = 0.1 /im, 0.5 /im and 3 /im in the case of Al we have Fi/Fi = 0.62, 0.89, 
0.98 from Eq. (Jg$ and 0.63, 0.88, 0.97 from [2J. For Cu and Au Eq. ® 
gives Fi/F[ 0) = 0.59, 0.85, 0.97 in agreement with the values: 0.55, 0.85, 0.97 
Bill. 



As was mentioned in the Introduction the computation of finite conductivity 
corrections to the Casimir force by the use of tabulated data was firstly 
performed in |*n| . It should be emphasized that our analytical results are in 
contradiction with JTDJ. By way of example, at a = 0.5 /im for Au and Cu one 
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can find in [|T9| F p /F^> = 0.657 and 0.837 correspondingly whereas according 
to our results F p /F^ = 0.81 for both metals. At the same distance and 
metals for a lens above a plate Fi/F^ = 0.719 and 0.874 whereas from 



Eq. ( gg| ) one gets Fi/Fi = 0.85. Our results, however, are in good agreement 
with the alternative computations of j2T| supporting the conclusion of [^T 



that the manipulation of optical data in [JT^] is invalid. 

It might be well to compare also the exact third and fourth order conductivity 
corrections obtained above with the approximate ones obtained by the use 
of interpolation formula [0, |TUfl . To take one example, for the force between 
a lens and a plate the coefficients near the third and fourth order corrections 
in interpolation formula are —50.67 and +177.33 (compare with —43.57 and 
+104.13 from (|2D|)). For the smallest separations a = 120 nm in experiment 



and 5o/a « 0.13 for Al this leads to the 0.5% difference only in the results 
obtained by the interpolation formula j7j and by 



5. Conclusions and discussion 

In the above the third and the fourth order corrections to the Casimir 
force due to finite conductivity of the metal were calculated analytically in 
configurations of two plane parallel plates and a spherical lens (or a sphere) 
above a plate. The Casimir forces (|24p, (|Zgf) are in excellent agreement with 
computations of |2Tj] based on the tabulated data for the complex dielectric 



permittivity for all distances larger than the effective plasma wavelength of 
the test body metal. What this means is that the results ( |Z4]) , ( |2D| ) can be 
reliably used even for the distances a less than the characteristic absorption 
wavelength Ao if X p < Ao (this is a case, e.g., for Au and Cw, which are 
characterized by Ao ~ 500 nm or for Cr with X p « 314 nm, Ao ~ 600 nm j30j). 

To obtain the higher order conductivity corrections we have used the 
plasma model representing dielectric permittivity by the Eq. (|I5|). This model 
does not take into account relaxation processes. However, the relaxation 
parameter is much smaller than the plasma frequency. As was shown in 
|2"Tf , relaxation could play some role only for large distances between plates 



a ^> Xp and even there the variation of the corrections to the Casimir force 
due to it is smaller than 2%. If to take into account that for so large distances 
the corrections themselves decrease very quickly it becomes evident that the 
influence of relaxation can be neglected. 

In conclusion we would like to stress that both the results ( |24p and 
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are of the same accuracy in spite of the fact that the PFT was used in 
(jT5p to obtain (p9|). The thing is that this theorem is equivalent to the 
addition method of calculation of the Casimir force which leads to the error no 
larger than 10 -2 % for small deviations from plane parallel geometry [|T2| , |3T 



What this means is that for a sphere or spherical lens of large curvature 
radius R^$> a the additional error introduced by the use of PFT is negligible. 
Therefore it is possible to conclude that formulas like ([Z3|), flZgp and the 
analogical expressions for the other corrections to the Casimir force can be 
reliably used for confronting theory and experiment at a level of 1% accuracy. 
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Figure captions 

Fig.l. Correction factors to the Casimir force F p /F^ in configuration 
of two plane parallel plates for Al (a) and Cu or Au (b) bodies 
in dependence of distance measured in /im. Solid lines represent 
the results of computations [21], short- and long-dashed lines are 
obtained by the Eq. flZ3Q up to the fourth and the second power 
respectively. 

Fig. 2. Correction factors to the Casimir force Fi/Fjy in configuration of 
a lens (sphere) above a plate for Al (a) and Cu or Au (b) bodies 
in dependence of distance measured in /im. Solid lines represent 
the results of computations [21], short- and long-dashed lines are 
obtained by the Eq. (|29|) up to the fourth and the second power 
respectively. 
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Figure l,a: 
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Figure l,b: 
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Figure 2, a: 
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Figure 2,b: 



